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We present results on the self-energy correction to the energy levels of hydrogen and hydrogenlike 
ions. The self energy represents the largest QED correction to the relativistic (Dirac-Coulomb) 
' energy of a bound electron. We focus on the perturbation expansion of the self energy of non-S 

, states, and provide estimates of the so-called Ago perturbative coefficient, which can be considered 

as a relativistic Bethe logarithm. Precise values of Aao are given for many P, D, F and G states, 
f~| ' while estimates are given for other electronic states. These results can be used in high-precision 

, spectroscopy experiments in hydrogen and hydrogenlike ions. They yield the best available estimate 

of the self-energy correction of many atomic states. 
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I. INTRODUCTION 

The recent dramatic progress in high-precision spec- 
troscopy (see, e.g., |1|) has motivated the calculation of 
numerous contributions to the energy levels of hydrogen 
and hydrogenlike systems. Such spectroscopic experi- 
ments test our understanding of atomic levels, and pro- 
vide precise determinations of fundamental constants ; 
this requires accurate predictions of atomic energies, and, 
in particular, the calculation of corrections due to Quan- 
tum Electrodynamics (QED), the quantum field theory 
of electromagnetic interactions. The largest correction 
to the relativistic (Dirac) energy levels of hydrogen and 
hydrogenlike ions is provided by the so-called self-energy 
contribution of QED. The self energy is a process which 
modifies the relativistic (Dirac) energy of an electron, 
and which can be depicted by the following Feynman di- 
agram. 



where the double line denotes the electron (bound to the 
nucleus), and where the wavy line represents the photon 
emitted and reabsorbed by the electron. The self-energy 
correction to energy levels in hydrogen and hydrogenlike 
ions can be expressed as an expansion in Za and \n(Za) 
(see, e.g., Q) — Z is the nuclear charge number of the 
nucleus of the hydrogenlike ion under consideration, and 
a is the fine-structure constant. Analytic calculations of 
the (one-loop) self energy in bound systems have a long 
history, starting from Bethe's seminal paper and have 
since extended over more than five decades. 

The purpose of this paper is to provide good approx- 
imate values of the self-energy correction to the energy 
levels of hydrogen and hydrogenlike ions, for any P state, 
and any state with a higher angular momentum. Only 
a part of the perturbation expansion of the self energy 
of these states is known analytically. The first two non- 
analytically-known contributions to this expansion are 
the Bethe logarithm lnA:o("'0 ^^'^ the so-called Aeo{nlj) 



coefficient of the self energy, which can be characterized 
as a relativistic Bethe logarithm [see Sec. ^1 and in par- 
ticular Eqs. Q and ©]. Here, nlj is the standard 
spectroscopic notation for an atomic state. This paper 
thus contains formulas for estimating both of these im- 
portant quantities (see Sec.lVland I VII) . 

Very precise numerical values of the Bethe logarithm 
lnfco("-0 have been obtained (see, e.g., Refs. Q and 
and numerical convergence acceleration techniques 
can yield very precise values of this quantity for any 
atomic state nl. The estimate (|37|) that we obtained as 
a by-product in Sec. IVII should be useful to experiments 
that use levels for which no published values of the Bethe 
logarithm exist (see, e.g., Ref. 8]). 

Many new values of the relativistic Bethe loga- 
rithm AQoinlj) have recently been published @. Other 
values have been obtained previously for some S 0, 0, 
and P states 0, Q . This paper contains two ad- 
ditional values [^60(5^5/2) and Aqo{5F'j/2)], as well as 
details on the procedure that we used in obtaining the 
values of Ago in Ref. ^ and in Table UTIl (see Sec. El- 

The results of Sec. liVHVil provide an improvement over 
the available approximations of the bound-electron self 
energy, over a large range of nuclear charge numbers Z. 
In particular, they yield the best available estimates for 
the self-energy correction in hydrogen, for all the states 
for which no exact (non-perturbative) value of the self 
energy has yet beeiijDublishcd (i.e., all levels, except n — 
1 and n = 2 levels [3il5|). 

It is important to know accurately the energy (and in 
particular the self energy) of higher angular momentum 
states, because they are used in high-p recision spectro- 
scopic measurements [wL UJl flSl Il9l l20l l2l| . States with 
very-high angular orbital quantum numbers I ~ 30 have 
been recently used in such experiments 's'l . Further moti- 
vation for the present study results from the need to accu- 
rately compare the two approaches that have been used 
for the theoretical study of QED shifts, so as to check 
their consistency: (i) the analytic expansion in the pa- 
rameter Za, mostly used for \ow-Z systems, and (ii) the 
numerical approach, which avoids the Za expansion and 
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has been used predominantly for the theoretical descrip- 
tion of high-Z hydrogenlike ions |22j| . 

Recently, the most accurate methods implementing a 
non-perturbative calculation of the self energy [isl 
0, I25l .26] have been extended by analytic results |27j . 
Taken together, they provide access to the self-energy 
shift of electrons of total angular momentum j > 3/2. 
This has allowed us to obtain numerical values of the self 
energy, and to use them in checks of the Ago coefficients 
presented in Tables HHTvl f see Sec. lVlUl . 

Moreover, general progress in theoretical calculations 
of atomic energy levels has been achieved by means of nu- 
merical algorithms 0, HE Hll that lead to an accelerated 
convergence of the angular-momentum series expansion 
of the bound-electron relativistic Green function. Such 
algorithms are also useful for performing the series sum- 
mations that we had to do in order to obtain the values 
of Aqo presented here (see Sec. II V|) . 

Notation and conventions are defined in Sec. ^ The 
mathematical method used for the semi-analytic calcula- 
tions of Ago in Ref. is discussed in Sec. IIIII Details 
on these calculations are presented in Sec. IIVI Formu- 
las for the relativistic Bethe logarithm AQo{nlj) of P and 
D states are presented in Sec.0 Estimates of the Bethe 
logarithm hiko{nl) and of AQQ{nlj) as a function of the 
orbital quantum number I are reported in Sec. IVll We 
have performed additional checks of the values of Aqq in 
Tables Hl lIVI as described in Sec. I VIII we also show in 
that section that for the states considered here, the in- 
clusion of in the (truncated) perturbation expansion 
of the electron self energy [Eq. below] does indeed 
improve the self energy estimates. A summary of the 
paper is given in Sec. IVIlH The fitting method that we 
used in obtaining asymptotic behaviors of Inko^nl) and 
of Aeo{nlj) is described in the Appendix. 



II. NOTATION AND CONVENTIONS 

In this section, we define the notation and conventions 
used in this paper. We write the (real part of the) one- 
loop self-energy shift of an electron in the level n with 
orbital angular momentum I and total angular momen- 
tum j as 

AEsE^-^-^F{nl„Za)mc\ (1) 

where F{nlj,Za) is a dimensionless quantity. We use 
natural units, in which fi = c = m = 1 and = Aira (m 
is the electron mass). It is customary in the literature to 
suppress the dependence of F on the quantum numbers 
n, j and / and write F{Za) for F{nlj, Za). 

The quantum numbers I and j can be combined into 
the Dirac angular quantum number k. As a function of 
j and K is given by 

n = 2{l~-j){3 + l/2), (2a) 



i.e., 

n = -{j + l/2) for i = / + l/2, (2b) 

and 

K = (j + l/2) for j = /-l/2. (2c) 

The quantum numbers j and I can be derived from k 
according to 

/ = |k+ 1/2| - 1/2 (3) 

and 

3 = kl - 1/2, (4) 

i.e., K specifies uniquely both j and I. The semi-analytic 
expansion of F(nlj, Za) about Za = for a general 
atomic state with quantum numbers n, / and j gives rise 
to the expression 

F{nlj,Za) = Aii{nlj)\n[{Za)-'^] 

+AiQ{nlj) + {Za) A5o{nlj) 
+ {Za)^ {A62{nlj) \n^[{Za)-^] 
+Aei{nlj) ln[(Za)-2] 
+GsE{nlj,Za)}. (5) 

This expansion is semi-analytic, i.e., it involves powers 
of Za and of ln[(Za)^^]. Terms added to the leading 
order in Za are commonly referred to as the binding 
corrections. The A coefficients have two indices, the first 
of which denotes the power of Za [including those powers 
contained in Eq. (Q], while the second index denotes the 
power of the logarithm \n{Za)^^. 

The limit as Za — *■ of GsEinlj, Za) is known to be 
finite and is referred to as the Aqq coefficient, i.e., 

Aeoinlj) = lim GsE{nlj,Za). (6) 

Za — >0 

Historically, the evaluation of the coefficient Aqq has 
been highly problematic. Due to the large number of 
terms that contribute at relative order (Za)'^ in (jSJ and 
problems concerning the separation of terms that con- 
tribute to a specific order in the Za expansion, eval- 
uations are plagued with severe calculational as well 
as conceptual difficulties. For example, the evaluation 
of AgnflSi /t) has drawn a lot of attention for a long 
time 0, 0, 113 113, 1^ • In general, the complexity of the 
calculation increases with increasing principal quantum 
number n. 

For many states, some of the coefficients in Q vanish. 
Notably, this is the case for P states and for states with 
higher angular momenta, as a consequence of their be- 
havior at the nucleus, which is less singular than that of 
S states (specifically, we have AQ2{'nlj) = A5o{nlj) = 
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A4i{nlj) = for Z 7^ — see Refs. 0, and refer- 
ences therein). The fact that the logarithmic coefficient 
Ari{nlj) contained in GsF.fn/,. Zq) in (0) vanishes for 
I ^ has been pointed out in 
pected that Aikinlj) 
thus have 



Za) i 

'32]; 



it is therefore ex- 
for k > 1. For nonzero I, we 



F{nlj,Za) = A4o{nlj) + (Zaf [A6i(n/j) ln(Za)"2 
+Aeo{nlj)]+0[iZaf] (1^0). (7) 



For the comparison to experimental data, it is useful to 
note that the terms in (0) and (0 acquire reduced-mass 
corrections according to Eqs. (2.5a) and (2.5b) of [33l |. 

The general formula for A40 for a non-S state reads 
(see, e.g., iilH) 



Aiainlj) 



1 



2k {21 + 1) 3 



--lnA:o(^0. (8) 



where the Bethe logarithm Infco(nZ) is an inherently non- 
relativistic quantity, whose expression reads p?4i § 19] 



In ko{nl) 



X 



2{Za)^m 



(9) 
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tained in the quoted equation. We have 



A(ix{nlj) 



(l-<5,,o) (3-^^^ 



n (21 + m) 

m— — 1 



(10) 



1 1 



601 77 



240 60 n2 

+7 In 2 + 3 (7 - In n + *(« + 1)) 



Herej_4' denotes the logarithmic derivative of the F func- 
tion In § 6.3], and 7 is Euler's constant [H § 6.1.3]. We 
may infer immediately 



^6l("-Pl/2) 
-461 ("•^'3/2) 

Aiii{nlj) 



1 / 29 

33- 



2 7 



45 
32 ( 3 



/(/ + 1) 

„2 



3 n {21 + m) 

m— — 1 



(lla) 
(lib) 

(/ > 2). (11c) 



For a given orbital angular momentum the coefficient 
Agi approaches a constant as n — )■ 00. Equation (|llc|) 
implies that A^\ is spin-independent for I > 2, i.e., for D, 
F, G, . . . states. Therefore, Aqi does not contribute to 
the fine structure of these states. 



Here, is the nonrelativistic Coulomb Hamiltonian 
p2/(2m) — {Za)/r, are the components of the mo- 
mentum operator {i is summed over from 1 to 3), and 
the ket 10) represents the Schrodinger wavefunction of 
a state with quantum numbers {n,l), with associated 
bound-state energy En — — {Z a)^ m / {2 n^) . The Bethe 
logarithm is spin-independent and therefore independent 
of the total angular momentum j for a given orbital an- 
gular momentum /; it can be written as a function of n 
and I alone [factors of Z cancel out in Eq. (O , so that the 
Bethe logarithm does not depend on Z] . For the atomic 
levels under investigation here, the Bethe logarithm has 
been evaluated in Refs. jESIHlIM EElEili3,lilli2 
(the results exhibit varying accuracies). Because Aqq in- 
volves relativistic corrections to the coefficient A^q , which 
in turn contains the Bethe logarithm, it is natural to refer 
to as a "relativistic Bethe logarithm." 

A general analytic result for the logarithmic correction 
Agi as a function of the bound state quantum numbers 
n, I and j can be inferred from Eq. (4.4a) of 3, 29] upon 
subtraction of the vacuum-polarization contribution con- 



III. THE e METHOD 

In this section, we illustrate the usefulness of the so- 
called e method [ill Ha . [lj| in bound-state calculations 
of QED corrections. It is known that relativistic correc- 
tions to the wavefunction and higher-order terms in the 
expansion of the bound-electron propagator in powers of 
Coulomb vertices generate QED corrections of higher or- 
der in Za (see, e.g., Ref. |4j and references therein); 
these terms manifest themselves in Eq. (O in the form 
of the function GsE{nlj, Za), which summarizes these ef- 
fects at the order of a {Zotf" m — see Eqs. and 0. It 
is also well known that for very soft virtual photons, the 
potential expansion fails and generates an infrared diver- 
gence, which is cut off by the atomic momentum scale, 
Za (see, e.g., Ref. |44,] and references therein). This cut- 
off for the infrared divergence is one of the mechanisms 
that lead to the logarithmic terms in Eq. jSJ. 

The e method is used for the separation of the two 
different energy scales for virtual photons: the nonrel- 
ativistic domain, in which the virtual photon assumes 
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values of the order of the atomic binding energy, and the 
relativistic domain, in which the virtual photon assumes 
values of the order of the electron rest mass. We con- 
sider here a model problem with one "virtual photon," 
that involves the separation of the function being inte- 
grated into a high- and a low-energy contribution. This 
requires the temporary introduction of a parameter e; the 
dependence on e will cancel at the end of calculation [see 
Eq. H22|) below] when the high- and the low-energy parts 
are added together. We have, 



nonrelativistic domain ^ e ^ electron rest mass. 



[Zaf 



m <^ e ^ m. 



(12) 



The high-energy part is associated with photon energies 
u; > e, and the low-energy part is associated with photon 
energies to < e. 

In order to illustrate the principles behind the 
e method, we discuss a simple, one-dimensional exam- 
ple: the evaluation of 



J{Za) 



1 



: dw, 



(13) 



where the integration variable uj may be interpreted as 
the "energy" of a "virtual photon." The integral J can be 
explicitly calculated, so that the perturbation expansion 
can be checked: 



2 {Zaf In 



1 



{Za) 



[^l-{Zar + \) 



J{Za) = -- 

For \Za\ < 1, this formula is uniquely defined; for other 
values of Za, the analytic continuations of the logarithm 
and of the square-root have to be performed consistently 
with the original definition H13|) . 

Within the e method, we start by dividing the calcu- 
lation of J{Za) into a high-energy part Jn^Za, e) and a 
low-energy part Ji^{Za, e), each of which depends on an 
additional parameter e [that satisfies (|12|) ]. The sum of 
the high- and low-energy contributions, which is 



J{Za) = Jii{Za, e) + Ji,{Za, e) 



(15) 



does not depend on e. Thus, the dependence on e should 
vanish entirely when we add the high- and low-energy 
contributions. We may therefore expand both contribu- 
tions Jh and Jl first in Za, then in e, and then add them 
up at the end of the calculation in order to obtain the 
semi-analytic expansion of J[Za) in powers of Za and 
In(Za). 

Let us first discuss the "high-energy part" of the cal- 
culation. It is given by the expression 



J-a{Za,e) 



1 



:dLU, (16) 



{Zaf+uj 

where it is important to note in particular the lower in 
tegration limit, e. For lo > e, we may expand 

(Za)^ - UJ 



{Zaf 



-1 



2 {Zaf 



UJ 



+ 0[{Zay 



(17) 



[see Eq. (|12() with m = 1]. Each corresponding term 
of (|16(l can be integrated, with result 



JuiZa^e) : 

TT 

2 



2{Za) 



In 



(18) 
■0[{Za)% 



where the ellipsis represent terms that vanish as e 0. 
It is sufficient to only include terms that don't vanish 
as e — > 0, to each order in Za, because the sum J in 
Eq. IjlSfl does not depend on e. Moreover, this makes the 
calculation more manageable. The full cancellation of 
the dependence on Ine will be explicit after we evaluate 
the "low-energy part." 

The contribution of the low-energy part (0 < w < e) 
reads 



JhiZa, e) 



{Zaf - UJ 1 
{Zaf+uj y/T^ 



■.duj, (19) 



where the upper limit of integration depends on e. For 
UJ < e, we use an expansion that avoids the infrared di- 
vergences that we encountered in Eq. H17|l : 



(20) 



which leads to a Za expansion of the low-energy part. 
We obtain for Jl: 



Ji,{Za,e) = 
{...) + 2 {Zaf 



In 



{Zaf 



(21) 

0[{ZafW{Za)], 



where the ellipsis again represents terms that vanish 
as e — > 0, and where j is some integer. 

When the high-energy part (|18|l and the low-energy 
part (|21() are added, the logarithmic divergences in e can- 
cel, as it should, and we have 

J{Za) — J-ii{Za,e) + J-L{Za,e) 

+ 2 {Zaf (ln[(Za)-2] +ln2) 



TT 

2 



-0[{ZafW{Za)] 



(22) 



(for some j), which is consistent with H14|) . We note the 
analogy of the above expression with the leading-order 
terms of the Za expansion of the function F{nlj, Za) 
given in Eq. iQ for I ^ (terms associated to the co- 
efficients A/^Q, Agi, and A^q). In an actual Lamb shift 
calculation, the simplifications observed between terms 
containing e are crucial 13, x\ . 

In this model example, the epsilon method allowed us 
to obtain H22() with minimal effort. For comparison, the 
reader may consider App. A of |45j , which illustrates the 
cancellation of e in higher orders of the Za expansion, 
using a different example. 
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IV. CALCULATION OF SELF-ENERGY 
COEFFICIENTS 

This section, along with the previous one, gives detail 
on the methods we used in order to obtain the values 
of the Ago coefficient in Tables Hl jlVI (see also Ref. 0). 
The purpose of our calculations is to provide data for 
the self-energy coefficients up to and including the rela- 
tive order {Zotf' [see Eq. 0]; for the states of interest 
here (non-S states) this corresponds to the coefficients 
A40, Agi and Ago- Equation (jSJ is the well-known gen- 
eral formula for the coefficient A40. The coefficient 
can be found in Eq. (|f 0(1 . with special cases treated in 
Eqs. (|lla|l - ()llcf) . The remaining nonlogarithmic term 
^60 is by far the most difficult to evaluate, and the first 
results for any with orbital angular momentum quantum 
number I > 2 were recently obtained in Ref. by using 
the methods described in this section. 

As explained in detail in 0, 0, 0| , the calculation of 
the one-loop self energy falls naturally into a high- and a 
low-energy part {Fh and Fl, respectively). In Sec. IIIII 
we illustrated this procedure, and the introduction of 
the scale-separation parameter e for the photon energy. 
According to Eqs. (39)-(43)], the contributions to 
the low-energy part can be separated naturally into the 
nonrelativistic dipole and the nonrelativistic quadrupole 
parts, and into relativistic corrections to the current, to 
the Hamiltonian, to the binding energy and to the wave- 
function of the bound state. We follow here the approach 
outlined in Refs. with some modifications. 

One main difference as compared to the evaluation 
scheme described previously concerns the nonrelativis- 
tic quadrupole (nq) part. It is given by a specific matrix 
element (see the definition of Pnq in Ref. |13l Eq. (39)]), 
which has to be evaluated on the atomic state and aver- 
aged over the angles of the photon wave vectors: 



reads 



An 



47r 6m 



i ik ■ r 

p e 



P 



1 



Hs~{E~ uo) 
1 



p' e 



-ih • r 



■ p' 



Hs-{E- uj) 
where the transverse 5 function is given by 



(23) 



The dipole interaction obtained by the replacement 

exp(i k ■ r) ^ 1 

is subtracted; it leads to a lower-order contribution. The 
next term in the Taylor expansion of the exponential 



Att 6m 



(24) 



(fc • r) 



1 



Hs-{E- uo) 
1 



pp {k ■ r) 



P 



pp (fc • r) 



This representation makes an evaluation in coordinate 
space possible. However, an evaluation of this expression 
leads to a rather involved angular momentum algebra. 
Specifically, we employ a well-known angular momentum 
decomposition of the coordinate-space hydrogen Green 
function 

G{ri,r2,E - uj) = ^ 5// (ri, rj, J/) Yi',™ (ri) (r2) , 

I' ,m 

(25) 

with E - u} ^ -~o? m/(2;/2), and jig 



4m /2riY /2r2 
av \ av \ av 



oc t21' + 1 f2ri\ r2l' + l /■2r2\ 

k=a 



(fc + l)2i'+i (I' + l + k-u) 



(26) 



where a = l/{Zam), (fc)c is the Pochhammer symbol, 
and L denotes associated Laguerre polynomials [43 . For 
a reference state |0) of orbital angular momentum I, we 
obtain in H24() nonzero contributions from Green-function 
components (jSH)) with I' ^ I - 2,1 ~ 1,1,1 + 1,1 + 2. They 
can be obtained by straightforward, but tedious applica- 
tion of angular momentum algebra (see, e^^ E3)- 

As in previous calculations (see also [Ij, Eqs. (18) 
and (19)] and 13, Eqs. (55)-(58)]), we obtain for the 
high-energy part of all atomic states the general struc- 
ture 



Fninlj, Za) 



1 



2k (2/ + 1) 



(27) 



+ ... 



Aei ln(2e) + 0{e) 



where /C is a constant, and where the ellipsis denotes 
higher-order terms [in Za and In(Za)]. As observed 
in Sec. IIIII we may suppress terms that vanish in the 
limit e ^ [terms of the form 0{e) in the (ZQ;)^-term in 
Eq. 1(27(1 above] . These terms cancel when the high- and 
low-energy parts are added. 

Together with the constant term — Agi In 2, the con- 
stant /C contributes to Aqq. C is the coefficient of the 
1/e divergence; the term —C/e cancels when the high- 
and low-energy parts are added. Both /C and C are state 
dependent and vary with n,j,l. As in 13, Eqs. (56) 
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and (57)] and "li, Eqs. (89)-(92)], the low-energy part, 
for all states under investigation, has the general struc- 
ture 



FL{nlj,Za) 



(Zaf 



; In kf){nl) 
e 



(28) 



Aei In 



0{e) 



where Infco(nZ) is the Bethe logarithm [see Eq. Q], and 
where the ellipsis denotes higher-order terms. The can- 
cellation of the divergence in e between (|27f) and H28|) 
is obvious. The constant C, which is state-dependent (a 
function of n, j, Z), represents the low-energy contribution 
to Aqo and can be interpreted as the relativistic general- 
ization of the Bethe logarithm. In terms of the general 
expressions (j7f|) and (|2Hl, ^eo is therefore given by 



IC-Aei \n2 + C. 



(29) 



Our improved results for Aqq coefficients rely essen- 
tially on a more general code for the analytic calcula- 
tions, written in the computer- algebra package Mathe- 
matica ,42, ^HJ , which enables the corrections to be eval- 
uated semi-automatically. Intermediate expressions with 
some 200,000 terms are encountered, and the complexity 
of the calculations sharply increases with the principal 
quantum number n, and, as far as the complexity of the 
angular momentum algebra is concerned, with the orbital 
angular quantum number of the bound electron. 

Of crucial importance was the development of conver- 
gence acceleration methods which were used extensively 
for the evaluation of remaining one-dimensional integrals 
which could not be done analytically. These integrals are 
analogous to expressions encountered in previous work 
(see [il Eqs. (36), (47) and (48)] and 14, Eqs. (80)~ 
(84)]). The numerically evaluated contributions involve 
slowly convergent hypergeometric series, and — in more 
extreme cases — infinite series over partial derivatives of 
hypergeometric functions, and generalizations of Lerch's 
<I> transcendent |^ [H^l ■ As a result of the summation 
over V in l|25|l . after performing radial integrals, two spe- 
cific hypergeometric functions enter naturally into the ex- 
pressions for the bound-state matrix elements that char- 
acterize the one-loop correction (see, e.g., jQ, Eqs (80) 
and (81)]). One of these functions is given by 



1-t^ 
1 + 



(30) 



where the integration variable t is in the range 0-1, and 
n is the bound-state principal quantum number (2F1 de- 
notes the hypergeometric function — see, e.g., Chap. 15 
in Ref. j42i]). For t ~ 0, the power series expansion of $1 
is slowly convergent. 



00 f 



2k 



k=0 



nt 



(31) 



TABLE I: Self-energy coefficient Aao © for P states [see 
Eq. O]. The quoted error is due to numerical integration. 
As in previous calculations (see Refs. [l^ . [T3 '). certain re- 
maining one-dimensional integrals involving (partial deriva- 
tives of) hypergeometric functions could only be evaluated 
numerically. 



Pl/2 (^=1) 



-0.998 904 402(1) 
-1.148 189 956(1) 
-1.195 688 142(1) 
-1.216 224 512(1) 
-1.226 702 391(1) 
-1.232 715 957(1) 



3/2 



(^ = ^2) 



-0.503 373 465(1) 
-0.597 569 388(1) 
-0.630 945 795(1) 
-0.647 013 508(1) 
-0.656 154 893(1) 
-0.662 027 568(1) 



The series is nonalternating. In order to accelerate the 
convergence in the range t £ (0,0.05), we employ the 
combined nonlinear-condensation transformation 7, '28| . 
The other hypergeometric function that occurs naturally 
in our calculations is 



$2(?^, t) = 2F1 1, ~nt, 1 - nt, 



1 - t 

T+t 



(32) 



For < t < 0.05, we accelerate the convergence of the 
alternating power series 



k=0 



nt — 



(33) 



via the S transformation jsj, Eq. (8.4-4)]. The conver- 
gence acceleration leads to a much more reliable eval- 
uation of the remaining numerical integrals which con- 
tribute to Aqq (but cannot be expressed in closed ana- 
lytic form). As a by-product of our investigations, we 
obtained through this (independent) method Bethe log- 
arithms which are consistent with the precise results of 
Ref. 5I]. Here, we restrict the accuracy to 24 figures and 
give results for P states: 

lnfco(2P) = 

-0.030 016 708 630 212 902 443 676(1), 

lnfco(3P) = 

-0.038 190 229 385 312 447 701 163(1), 

lnfco(4P)= (34) 
-0.041 954 894 598 085 548 671 037(1), 

lnfco(5P) = 

-0.044 034 695 591 877 795 070 318(1). 

These results, which test the numerical methods that we 
employed, are in agreement with other recent calcula- 
tions IE la El El 

The main results of this paper concerning the Aqq co- 
efficients are given in Tables U HlVI with an absolute pre- 
cision of 10~^. In addition, we give explicit expressions 
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TABLE 11: Aeo coefficients © for D states. 



n 

X 

4 
5 
6 
7 

8 



D3/2 = 2) 



D5/2 



- -3) 



0.005 551 573(1) 
0.005 585 985(1) 
0.006 152 175(1) 
0.006 749 745(1) 
0.007 277 403(1) 
0.007 723 850(1) 



0.027 609 989(1) 
0.031 411 862(1) 
0.033 077 570(1) 
0.033 908 493(1) 
0.034 355 926(1) 
0.034 607 492(1) 



TABLE IIL Aeo coefficients © for F states. 



n 

T 

5 



F5/2 = 3) 



0.002 326 988(1) 
0.002 403 158(1) 



7/2 



{k = -4) 



0.007 074 961(1) 
0.008 087 020(1) 



for the low- and high-energy parts of the self energy, for 
the states with n — 5 under investigation [see Eqs. H27(l 
and (HHI) and Table E). They may be helpful in an in- 
dependent verification of our calculations. Note that the 
G7/2 and Gg/2 states involve the most problematic an- 
gular momentum algebra of all atomic states considered 
here. 

For some P states (see Tabled, the values of Ago re- 
ported here are four orders of magnitude more accurate 
than previous results 0, 0| , due to the improved nu- 
merical algorithms. For the 3Pi/2 states, the numerical 
value for the Ago coefficients of Table |2 differs from the 
previously reported result |14| by more than the numer- 
ical uncertainty quoted in Ref. whereas they are in 
agreement with previous results[l3L [T3| in the case of 
2P1/2 and 4Pi/2 states. The discrepancy for A6o(3Pi/2) 
is on the level of 5 x 10~^ in absolute units, which cor- 
responds to roughly 2 Hz (in frequency units) on the 
self-energy correction in atomic hydrogen. The com- 
putational error in Ref. [l^ was caused by numerical 
difhculties in one of the remaining one-dimensional in- 
tegrals involving the hypergeometric functions H3U|) and 
(|32|l . which could not be evaluated analytically. The nu- 
merical difficulties encountered in previous calculations 
due to slow convergence of the integrals are essentially 
removed by the convergence acceleration techniques. 

For some states, rather severe numerical cancellations 
are observed between the high- and low-energy contribu- 
tions to the self energy, as well as between the different 
contributions to the low-energy part. This intriguing ob- 
servation is documented in Tables IVTl and IVTll using the 



5G7/2 state as an example. Note that these numerical 
cancellations go beyond the required exact, analytic can- 
cellation of the divergent contributions which depend on 



TABLE V: According to Eqs. ^ and the high- and 

low-energy parts can be cast into a general form involving 
the terms C, K, and C. The coefficient Ago can be expressed 
in terms of K, Aq\ and C according to 1291 1. Here, we present 
analytic results for the terms C, Ae,\ and /C, and numerical 
results for C (for states with n = 5). The results for Am can 
be inferred from Eqs. llOl - lllcll . For i > 2, we observe that 
the Aei are spin-independent and that C = Agi. 

C, /C and C coefficients for states with n = 5 

state C Abi /C C 



rp 292 796 

Ori/2 1125 1125 



5P3/2 
5D3/2 
5D5/2 
5F5/2 
5F7/2 
5G7/2 
5G9/2 



292 
1125 



92 
7875 



92 
7875 



2 

1125 



2 

1125 



2 

4725 



2 

4725 



436 
1125 



92 
7875 



92 
7875 



2 

1125 



2 

1125 



2 

4725 



2 

4725 



20129 
67500 

199387 
540000 

35947 
3780000 

3097 



157500 

2657 
1102500 

774121 
211680000 

4397 
6048000 

269 



283500 



-1.023 
-0.747 
0.023 
0.021 
0.006 
0.005 
0.001 
0.001 



991 781(1) 
615 653(1) 
759 683(1) 
511 798(1) 
045 397(1) 
662 248(1) 
834 827(1) 
757 471(1) 



TABLE VI: As explained in Refs. [HITl , the low-energy con- 
tributions to AdQ naturally separate into the following terms: 
the nonrelativistic quadrupole part 



relativistic corrections to the current Fg- 




1 



Eq. (39)], the 
Eq. (40)], rel- 
ativistic corrections to the Hamiltonian Fsyi \^ Eq. (41} 
and relativistic corrections to the bound-state energy Fss 
Eq. (42)] and to the wavefunction Fs^, 13, Eq. (43)]. This 
classification suggests that it is natural to refer to the low- 
energy contribution £ as a relativistic Bethe logarithm. Ob- 
serve that the total contribution to A^sq of the low-energy part, 
which reads 0.001 834 827(1), is roughly five times smaller 
than the largest individual contribution (from Fsn), due to 
mutual cancellations. 

Contributions to the low-energy part (5G7/2 



Aeo-contribution due to F^,^ 
Ago-contribution due to Fgy 
ylao-contribution due to Fsh 
A60-contribution due to Fse 
A60-contribution due to Fs^, 
A^o (see entry for £ in Table IVt 



0.002 875 830 9(5) 
-0.001 083 109 4(5) 
-0.008 917 782 1(5) 
0.004 920 556 0(5) 
0.004 039 332 1(5) 
0.001 834 827(1) 



TABLE IV: ^60 coefficients © for G states. 



n 


G7/2 = 4) 


G9/2 {1^ = -5) 


5 


0.000 814 415(1) 


0.002 412 929(1) 



the scale-separation parameter e. 
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TABLE VII: For the 5G7/2 state, an additional numerical 
cancellation occurs when the finite contributions to Ago origi- 
nating from the low-energy part (see the ninth row of TablelVt 
and the high-energy part are added according to Eq. (I29I I. 
The high-energy contribution is A6o(^h) = JC — Aei In 2, and 
the low-energy contribution is Aeo{Fi,) = £. 





-0.001 020 413 




0.001 834 828(1) 


Aeo 


0.000 814 415(1) 



TABLE VIII: The asymptotic behavior of Aeo{nlj) as n — > 00 
can be described by an expansion in 1/n. The following table 
contains the first coefficients of such an expansion, as defined 
in Eq. I|35^ . The approximate values of Aeo{nlj) that can be 
directly deduced from this table and from Eq. <35bll are the 
best available values of Aeo for P and D states, except for the 
states that are represented in Tables HI and HTl These results 
are depicted in Fig. 



state 


ao 


ffli 


a2 


Pl/2 


-1.249(9) 


0.0(2) 


0.87(45) 


P3/2 


-0.69(2) 


0.15(5) 


0.25(25) 


D3/2 


0.011(1) 


-0.032(7) 


-0.05(9) 


D5/2 


0.034(2) 


0.025(5) 


-0.18(4) 



energy of the electron of an hydrogenlike ion can be esti- 
mated through Eqs. Q, ©, l(TT|) a nd with A3 de- 
fined with the values of Table IVIIII Hydrogen has been 
and will be the subject of extremely precise spectroscopy 
experiments, which now reach the level of 1 Hz of uncer- 
tainty in transition frequencies. The uncertainty in the 
self energy that comes from the uncertainties in the 
coefficients of Table IVIIll through ||7J) and ^ IS com- 
parable to the current experimental limit. In fact, the 
uncertainties in A3 in (|35bl) contribute to the self energy 
by less than ±2 Hz for P1/2 states with n > 7, less than 
±1.6 Hz for P3/2 states with n > 7, less than ±0.12 Hz 
for D5/2 states with n > 8, and less than ±0.12 Hz for 
D5/2 states with n > 8 (precise values of Ago for lower 
values of n can be found in Tables III and llT|) . 

Moreover, the coefficients of Table IVIIll can be useful 
to theoretical calculations. In fact, future values of Ago 
for P and D states can be checked against the estimates 
provided by A3 in Ij35bp — see also the curves of Fig. ^ 



VI. APPROXIMATIONS OF Ago AND OF THE 
BETHE LOGARITHM 



V. Aeo FOR HIGHER-n STATES 

This section contains approximate formulas that we 
have found for the Ago coefficients of P and D states, 
for principal quantum numbers n that go beyond those 
of Tables |2 and These tables contain enough values 
of Ago(nlj) for extrapolations to be made. We present 
the asymptotic behavior of Ago(nZj) as n ^ 00 as 

Aeoinlj) = Asin, Ij) + O ^-^^ , (35a) 

where 

A3in, I,) = ao(/,) + ^ + (35b) 

Such an asymptotic behavior is justified, for any non- 
S state, by its similarity to the functional form of the 
self-energy coefficient Agi in Eq. {Tj) — see Eq. Hll|) . The 
values that we obtained for the coefficients ai{lj) can be 
found in Table IVTm The fitting method is described in 
the Appendix. 

The approximation A3(n, Ij) to Ago(ri?j) is depicted in 
Fig. for P and D states. According to the graphs in 
this figure, the 0{l/n?) contribution in H35a|) is much 
smaller than the uncertainty in A3, which comes from 
the uncertainties in the coefficients of Table IVIIII 

The coefficients of Ij35b|l given in Table IVIIII can 
be useful to spectroscopy experiments that involve elec- 
tronic levels with principal quantum numbers that are 
higher than those of Tables and ^ In fact, the self 



In addition to studying the dependence of Ago(nZj) 
on n, as we did in the previous section for P and D states, 
it is interesting to analyze the behavior of Ago {nlj ) as a 
function of I, for j = Z — 1/2 and j = I + 1/2. We con- 
jecture that Ago(nZj), ior n = I + 1 and j — I ± 1/2, 
decreases as 

c(? - I) 

AeoinL) - ^ ,, ^ with fc > 3, (36) 

where we probably have k = 4 or k = 5 [c(l/2) and 
c(— 1/2) are two unspecified numbers]. The form (|36|l is 
motivated in this section. 

We have also studied the asymptotic behavior of the 
Bethe logarithm Infco(nZ), because this is a quantity sim- 
ilar to the "relativistic Bethe logarithm" Ago, and be- 
cause it yields a large contribution to the self energy 
[see Eqs. {Tjl and JSJ]. We show in this section that the 
Bethe logarithm ln/co(^Oi where h = I + 1, appears to 
behave asymptotically as This result differs from 

the Z~^/^ asymptotic behavior of \nko{nl) deduced from 
Eq. (B5) in O p. 845]. Extrapolations of the Bethe 
logarithm In ko{nl) as a function of n have been obtained 
through the method described in the Appendix, and used 
in Ref. ^ for S, P and D states (Z = to 2). 

We also postulate that the Bethe logarithm \nko{nl), 
where h = I + 1, can be expanded in powers of 
about I — 00. In order to find the first five coefficients 
of such an expansion, we used the fitting procedure de- 
scribed in the Appendix. The resulting approximation 
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P 112 D3/2 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 

1/n 1/n 



FIG. 1: These graphs show exact and approximate values of the self-energy coefficient Aqq — see Eq. O. Exact values are 
represented by dots and can be found in Tables d and lUl The two curves of each graph represent the upper and lower limits of 
the approximation to Aeo provided by A3 in Eqs. 13511 . by taking into account the uncertainties in the coefflcients of Table IVTm 
For levels in hydrogen with principal quantum number n > 10, the uncertainty in ^eo deduced from these curves contributes 
to the uncertainty in the electron self energy Q by less than 2 Hz. (The use of 1/n as the abscissa allows all large principal 
quantum numbers n to be represented in the graphs.) 



reads: 



IK 



-0.046 
-0.048 
-0.05 
-0.052 
-0.054 
-0.056 
-0.058 



xlnkoihl) (37) 

/ ^^r.oro,„^ 0.02478(4) 0.0387(8) 
- 0.056853(2) + j—^ + j^-^ 

-0.114(6) 0.16(2)\ 
^ P + ;4 J' 

where h = / + 1, and where the neglected contribu- 
tion is of order l^^. This approximation should be valid 
for ^ — > 00; nevertheless, it yields values of the Bethe log- 
arithm that are both precise (see Fig. |2Jl and compatible 
with all the values of \nko{nl) for I = 3, . . . , 19 (taken 
from Ref. @). For the / > 20 levels of hydrogen, the un- 
certainty in the result of approximation H37|l is negligible, 
when compared to the best experimental uncertainty in 
transition frequency measurements (about 1 Hz 0]). 

Moreover, we suggest that the orders of magnitude of 
the self-energy coefficient AQo{nlj) and of the Bethe log- 
arithm lnfco("-0 do '^ot depend on the principal quan- 
tum number n, i.e., the order of magnitude of a co- 
efficient AQolnlj) is given by the order of magnitude 
of AQo{nlj), where n = I + 1 (and similarly for the Bethe 
logarithm). For Aqq, this behavior is a generalization 
of what is observed for P, D, F and G states in Tables 
H HIVI For the Bethe logarithm, the fact that \nko{nl) 
and lnfco("-0 have the same order of magnitude can be 
observed for states with I < n < 20 hy inspecting the 



0.1 0.2 0.3 0.4 

1-' 

FIG. 2; Comparison between exact values of x Infco(ni) 
(dots) and the truncated asymptotic expansion of Eq. 13 711 
(zone between the two curves) — lnA:Q(n?) is the Bethe loga- 
rithm, and n = 1 + 1. The numerical values of the Bethe 
logarithms used in this graph 6'] are compatible with the val- 
ues deduced from Eq. H37^ , which are in the area between the 
two curves. The fact that the data points seem to converge 
toward a finite value (~ —0.057) as supports the con- 

jecture of a asymptotic behavior of the Bethe logarithm 
In ko{nl). 



results of Ref. @ . 

The expressions H36|l and l|37|l for the asymptotic be- 
havior of AQo{nlj) and lnfco(fi/), where n = I + 1, could 
thus be used for estimating the order of magnitude of 
the self energy JlJ— with the help of Eqs. 0, ©, EJ. 
Estimating the self energy correction can be useful in 
high-precision spectroscopy experiments with large-Z lev- 
els. Thus, for instance, a recent experiment j^] required 
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evaluating the self energies of circular {n = I + 1) states 
of orbital quantum number I ~ 30. On the theoretical 
side, future calculations of AeQ{nlj) and Inko^nl) can be 
checked against the asymptotic behaviors of Aeo{nlj) and 
Inko^nl) that are described above. 

Since the order of magnitude of Aeoinlj) does not 
appear to depend on n, it is natural to represent it 
(for fixed I and j) by the order of magnitude of either 
lim„^oo ^6o('^'j) — largest possible n — or Aeo{hlj) — 
where n = l + l is the smallest n possible for the angular 
momentum quantum number I. We chose the latter pos- 
sibility for two reasons. First, small-n values of AeQ{nlj) 
are available (see Tables Ul lIVI) . Second, future values of 
Aeoinlj) for higher angular quantum numbers / are likely 
to be obtained first for states where n = I + 1, which is 
the smallest n possible for a given angular momentum 
quantum number I. In particular, such states have sim- 
pler radial wavefunctions (the number of terms in the 
radial wavefunction of a state increases with n — I). And 
finally, circular states {n — l + l) are relevant to high- 
precision spectroscopy experiments (see, e.g., Ref. JJ), 
whereas n = oo states are unphysical. 

As mentioned above, we expect an asymptotic behav- 
ior of the form with k integer, for AQo{nlj) and for 
the Bethe logarithm Infco(nl). Such a functional form is 
motivated by the fact that all the Aik{nlj) coefficients of 
the self-energy function F in Eq. (|SJ) can be expanded in 
power series of 1/n and except maybe for the two co- 
efficients related to this section, Ago and A^q, where the 
latter is a function of the Bethe logarithm [see Eq. ©]. 
(We suppose that Aqq and A40 can also be expanded in 
such a series.) This can for instance be checked with the 
formulas for Aik{nlj) reviewed in Ref. 0, p. 468], with 
the help of Eq. ifTUI) for AQi{nlj), where 5'(n + 1) can be 
expanded in powers of l/{n + 1) § 6.3.18]. 

The behavior of the Bethe logarithm \nko{nl), 
where n = l + l, is suggested by Fig. [3 The points 
of this graph, which represent 



P X lnfco("Oj 



(38) 



appear to converge toward a limit (~ —0.057) as 0. 
We checked the l~^ behavior deduced from the study of 
Eq. by calculating the slope of a log-log plot of the 
Bethe logarithm \iiko{'nl) (with numerical values taken 
from Ref. Q). The result, shown in Fig.O indicates that 
the Bethe logarithm does indeed behave asymptotically 
as l~^; this coincides with the conclusion from Fig.|3 

It is possible to use the procedure depicted in Fig. ^ 
to estimate the integer exponent k of an asymptotic be- 
havior l^'' for the relativistic Bethe logarithm AQo{nlj), 
where fi = / + 1 and j ~ I ± 1/2. In fact, it is rea- 
sonable to use the Bethe logarithm lnfco(nO ^ guide 
for studying the relativistic Bethe logarithm Aqq. Thus, 
the procedure depicted in Fig. |31 was applied to the self- 
energy coefficient AQo{filj); we obtained the asymptotic 
behavior presented at the beginning of this section, and 
in particular in Eq. H36|) . The graphs supporting l|36l) 
are given in Fig. ^ for with states with j = I + 1/2, and 
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FIG. 3: Upper graph: log-log plot of the Bethe logarithm 
lnfco(n/), where n = I + 1. Lower graph: slope between two 
successive points of the log-log plot. The limit slope of —3 
as / — > 00 observed in the lower graph indicates that the 
Bethe logarithm lnfco(nO behaves asymptotically as . This 
confirms what is observed in Fig. |21 (The abscissa of the 
points in the lower graph is chosen so as to produce a graph 
from which the limit slope of the upper graph as Z — > cx) can 
be easily deduced.) 



in Fig. [SI for states with j = I — 1/2. Each of these 
graphs uses only three values of Aqq (D, F and G states); 
even though this is a relatively small number of values 
compared to the number of available values of the Bethe 
logarithm, the behavior of the first few data points in 
Fig.|21justifies using only a few small-Z values in order to 
predict the asymptotic behavior of AQQ{nlj) for / — > 00. 

The values of the Ago coefficient of S and P states were 
not used in obtaining Eq. (|36|l . because it is convenient 
to treat the orders of magnitude of the Aqq coefficient of 
these states separately from the orders of magnitude of 
higher-Z states; Fig.Elillustrates this point. We note that 
the self-energy coefficient Aqi also exhibits an exceptional 
behavior for S and P states (see, e.g., Eq. (4.4a) in '^). 
As an additional consequence, estimating the coefficient c 
of the asymptotic form of Aqq in Eq. 136(1 would require 
use of states with orbital angular momentum quantum 
number l>2 {D, F, etc.). 

The possible values of the exponent k in Eq. H36|l de- 
duced from both the graphs of Fig. 01 and of Fig. El are 
compatible with each other (fc > 3 with, probably, fc = 4 
or fc = 5). It is indeed expected that the asymptotic form 
of AQo{nlj) be the same for j = I + 1/2 and j = Z — 1/2, 
as can be seen from the numerical values for D, F and G 
states found in Tables ITlljlVI More precise estimates of 
the asymptotic exponent k in Eq. H36|) can be obtained 
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0.1 0.2 0.3 0.4 0.5 0.60.7 0.8 

FIG. 4: Upper graph: log-log plot of the self-energy coef- 
ficient Afio{nlj), where n — I + 1 and j = I + 1/2. Lower 
graph: slope between two successive points of the log-log plot 
(solid line) and extrapolation to Z — > oo (dashes). By analogy 
with the graphs similarly obtained for the Bethe logarithm 
in Fig. 13 we conclude that for j — I + 1/2, A^oinlj) behaves 
asymptotically as with fc > 3 and, probably, fc = 4 or 
fc = 5. The values of Aeo are taken from Tables HTlll VI 



through the procedure we used in Figs. 0] and |S1 as soon 
as additional values of Aqq {nlj ) , with n = I + 1 are avail- 
able. 

According to the results of this section, the "relativistic 
Bethe logarithm" Aqq {nlj ) decreases at least as fast (and 
probably one or two powers faster), as a function of I, 
than the Bethe logarithm lnfco(?^0• Such a behavior is 
also found in the (Dirac-Coulomb) energy of hydrogen 
and hydrogcnlike ions. Thus, the Dirac-Coulomb energy 
of an electron bound to a nucleus of charge number Z is 
(see, e.g., 0, p. 466]) 



En 



(Za 



i2 1 



{n - (5)2 



-1/2 



(39) 



where 



5 ={2 + 1/2) - v/0- + 1/2)2 -(Za)2. 



According to (|39|l . an electron in a circular state nlj with 
j = Z + 1/2 (and n = Z + 1) has an energy 



(40) 



In the Taylor expansion (in Za) of this energy, the 
asymptotic behavior of the coefficient of [Zol)^^ is given 
by (this conclusion also holds for circular state filj 
with 2=1 — 1/2). Thus, for circular states, successive 
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FIG. 5: Upper graph: log-log plot of the self-energy coef- 
ficient Afjo(nlj), where n = / -|- 1 and j ~ I — 1/2. Lower 
graph: slope between two successive points of the log-log plot 
(solid line) and extrapolation to Z — > oo (dashes). By analogy 
with the graphs similarly obtained for the Bethe logarithm 
in Fig. 121 we conclude that for j ~ I — 1/2, A^oinlj) behaves 
asymptotically as with fc > 3 and, probably, fc = 4 or 
k — 5. The values of Aao are taken from Tables UlUlVI 
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FIG. 6: This graph shows values of the self-energy coeffi- 
cients AQ(i{nlj), where n = Z + 1, as a function of the Dirac 
quantum number k, where k is defined in ||5J . The large value 
^6o(l'S'i/2) — —31 is not represented here. This plot shows 
that for S and P states {k = —2, —1 and 1), the Aqq coefficient 
exhibits an exceptional behavior; such an exceptional behav- 
ior is also found in the self-energy coefficient Agi in Eq. llOt , 
which is known analytically. 



relativistic corrections to the nonrelativistic energy of a 
bound electron fall off faster and faster with the orbital 
quantum number Z, with two additional powers of Z^^ for 
each order in [Zaf' . If this rule applies to the coefficients 
of the self-energy expansion Q , the asymptotic form of 
AQQinlj) as Z — > oo should be Z~^; in fact, the lower-order 
coefficient A^Qifilj) decreases as Z"^, as can be seen in 
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Eq. ((SJ. On the other hand, since Aeo{nlj) can be con- 
sidered as a relativistic correction to the Bethe logarithm, 
applying the above rule yields an asymptotic form in 
for AeQ{nlj), since the Bethe logarithm behaves as l^^, 
as described in this section. These observations are fully 
compatible with the graphs of Figs. 01 and El from which 
the asymptotic form H3t)|) of AQQ^hlj) was deduced (with 
an exponent k probably equal to 4 or 5). 
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VII. CHECKS OF THE Ar,a COEFFICIENTS 

We have checked our analytic results for ^go (cf. 
Tables Hl lIVII by an independent method: the ana- 
lytic results were compared to values deduced from 
non-perturbative, numerical calculations of the self en- 
ergy Q). We have used the numerical self-energy values 

of Refs. HI mm mill, as well as new values 591 , 

which extend the results of Ref. jl^l to smaller nuclear 
charge numbers Z (to Z between 10 and 25). In most 
cases, the checks that we detail below confirm the values 
of reported in Tables ID lIVI to a relative precision of 
about 15 %. The few exceptions are the following. For 
2P states, the numerical values of the self-energy confirm 
the results of Tabled to about 1 %. For states 
with n — 3, . . . , 8, the non-pcrturbative self-energy re- 
sults yield AQo{nD^/2) — 0.005(10), in agreement with 
the results of Table m And finally, we did not check 
^4.60(805/2) in Table UTI by using non-perturbative self- 
energy values because no such values are available for the 
8D5/2 state. However, as depicted in Fig. ^ the value of 
^60(805/2) reported here appears to fit well within the 
series oi A^Q[nD^i2) values for n = 3, . . . , 7 (see Table lTHl . 

The first check that we applied consisted of comparing 
the numerical, exact results for F to two of its successive 
approximations. The first approximation, F'''^\Za), in- 
cludes the two dominant and already-known coefficients 
AiQ ^ and Aqi (^Ul of expansion (TJ: 



F^^\Za) = A40 + {ZafAei \n{Zay 



(41) 



and the second approximation, F'^^\ includes in addition 
the next-order contribution reported in this paper: 

{Za) = A40 + {Zaf [^61 ln(Za)-2 + Ago] • (42) 

For a given electronic level nlj, one expects that 
for low Z, the curve of the higher-order approxi- 
mation F''^\Za) be closer to the curve of F{Za) 
than F^'^'>{Za). In order to check this, we plotted the 
quantity 



I{Za) = In 



F{Za) - F^^\Za) 



F{Za) - F(2)(Za) 



(43) 



which should go to — cxd as Z — > 0, as can be seen from 
Eq. (0). In (|43|l . the purpose of the logarithm is only to 
obtain more legible graphs; a value of / lower than zero 
indicates that including Aqq in the approximation of F 
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FIG. 7: Figure (a) shows exact and approximate values of 
the (scaled) self energy F of a 2P1/2 electron [see Eq. Q]. 
Exact values are given on the solid line. The two-coefficient 
approximation 14111 is represented by long dashes. The three- 
coefficient approximation 14211 uses the value of j46o(2Pi/2) 
that we provide in Table HI and is indicated by short dashes. 
Figure (b) displays the improvement provided by the inclu- 
sion of Aso in the self-energy approximation, as measured by 
the function I in Eq. 1431 : negative values of I indicate that 
including A(,o improves the approximation. 



improves the lower-order approximation. For the states 
of Tables lIV llVI graphs of H43fl are compatible with their 
expected behavior [I{Za) is negative for Z sufficiently 
close to zero, and is consistent with a —00 limit]. Fig- 
ures [71 and jSl show this behavior for two electronic states. 

Moreover, the improvement provided by the inclusion 
of Aqq in the approximation for F becomes greater as the 
total angular momentum j increases: for given n and Z, 
the improvement function H43|) decreases as j increases; 
this behavior can observed by comparing Figs. [71 and |S1 
Similarily, the range of Z for which approximation i^'^' 
is better than F^^'' increases with increasing j. In the 
worst of the cases considered here (j — 1/2), approxima- 
tion F^'^) is better than F^^^ up to Z ~ 25. As shown 
in Fig. IHl for a high-j level such as 5G7/2, the higher- 
order approximation F^^^ is better than F^^^ even up 
to Z = 110. 

The second check consisted in estimating A^q from the 
numerical values of the self energy (Q. For all the elec- 
tronic levels nlj studied here (except for 8D5/2), we have 
plotted the function GsE{nlj, Za) of (jSJ; this is made 
possible by the fact that all the coefficients of (jSJ are 
(analytically) known for any state 0, except for 
the Bethe logarithm, which has been numerically eval- 
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FIG. 8: These two figures represent respectively the same 
quantities as those found in Fig. Q but for the 5G7/2 level 
instead of the 2P1/2 level. The fact that curve (b) contains 
negative values of I [see Eq. 103] indicates that the three- 
order approximation l|42^ to the self energy Q is better than 
the two-order approximation Mil , at least over the range of 
nuclear charge numbers Z = 25-110. The three-order ap- 
proximation Il42t uses the value of A(,o(J)Gr /2) reported in 
Table HVI 



uated for many states, including the ones we consider 
here 0, H IH [13. As indicated in ®, the hmit of 
the remainder GsEinlj, Za) as Za — > is by defini- 
tion AQQ{nlj). We have estimated this limit both visually 
and by fitting GsEinlj, Za) with various choices of non- 
zero higher-order terms. A typical curve for Gs^iZa) is 
shown in Fig.|Hl The estimates of Ago obtained by these 
procedures confirm the independent analytic results of 
Tables HI Hvl to a typical accuracy of 10-20 %, with a few 
exceptions. Thus, for 2P levels, plotting Gse as in Fig. El 
allowed us to confirm the values of AQo{2Pj) in Table |3 
to a precision of about 1 %. This higher precision is ob- 
tained by using the self energies of 2P states obtained in 
Ref. 15] for values of Za close to zero {Z = 1, . . . , 5): 
such low-Z self energies are well-suited to an evaluation 
of Ago by the limit ©. Plotting Gse for D3/2 states lead 
to AQo{nD:i/2) = 0.005(10) for n = 3,..., 8, in agree- 
ment with Table ^ Finally, since no non-perturbative 
self-energy Q is available for 8D5/2 states, we were not 
able to independently obtain ^00(805/2) by using such 
values. 

As a by-product of our work with graphs of 
GsEinlj, Za), we estimate the self-energy remainder 
GsE{nlj,a) relevant to hydrogen {Z = 1) to be 0.030(5) 
for 3D5/2 and 4D5/2 states [see Eq. this is larger than 
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FIG. 9: Plot (solid line) of numerical values of the remain- 
der Gse(4D5/2, -^a) of the self energy ||KJ; the dashed line 
indicates the value of ^60(405/2) — 0.0314 reported in this 
paper (see Table HTl. By definition, the coefficient Aqo can be 
obtained as the limit JHJ of Gse as Za — > 0. This plot shows 
that the value of Aao extracted from numerical self-energies 
is consistent with the value obtained by the calculations pre- 
sented in this paper. We made identical observations for all 
the states of Tables HHIVl 
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FIG. 10: Plot of the function K in Eq. ||I7|| for the 5F 



and 5F5/2 states. 



7/2 



The limit of this function as Z must be 
zero if the coefficients Aqq of Tables II1 41 V I agree with exact, nu- 
merical values of the self energy. The curve displayed here in- 
dicates that the two values of Afs^6o(5-F') in Eq. 14611 obtained 
independently from Table UTTl and from non-perturbative self 
energies Q [via Eqs. Q, 0, Ullcjl . and (0^1 do not 
differ by more than about 3 %. 



the estimate of 0.00(1) given in Ref. H p. 468]. These two 
new values change the previous estimate of the self energy 
of 3D5/2 and 4D5/2 states through Eq. {TJ by a relatively 
large amount, compared to the current best experimen- 
tal uncertainty in transition frequencies (about 1 Hz TJ). 
Thus, a variation of 0.03 in GsE(3D5/2,a) in (|SJ corre- 
sponds to a variation of about 50 Hz in the self energy 
correction ^ of the 3D5/2 level in hydrogen. The same 
variation in Gse(4D5/2, a) induces a variation of about 
20 Hz in the self energy of the 4D5/2 level in hydrogen; on 
the other hand, this latter change is small compared to 
the uncertainty of the relevant measurements considered 
in Ref. 0. 

As a third and last check, we used the numerical, ex- 
act values of F in order to study the following difference 
between remainders Gse [see Eqs. (01 and 0]: 



AfsGsE("-^ Za) = Gse(»^'(+i/2, Za) - GsE("-^i- 



-1/2 



Za), 
(44) 
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where, by definition of Aqq 

lim AfsGsE(?^^, ^a) 

Za^O 

= AtsA6o(nZ), 



(45) 
(46) 



which denotes a quantity associated to tlic fine-structure. 
The numerical evaluation of this limit is interesting: for 
the states of Tables U HlVI the numerical results for F 
yield values of AfsAeoinl) that are more accurate than 
our numerical estimates of the two individual terms 
^6o("-'i+i/2) and AQo{nli_i/2)- Our analytic values for 
AfeAeo in Eq. (03 were checked by plotting 



KiZ) 



A{sAeo{nl) 



1, 



(47) 



where AfsGsE('T-', Za) was calculated from the numerical 
values of F [see Eq. (0 and the coefficients reproduced 
in Sec.m, and where the value of AfsAQo{nl) in Eq. ll^ 
was deduced from the analytic results of Tables U HlVI 
If the numerical and analytic estimates of A{sAQo{nl) do 
agree, the function (|47|l goes to zero as Z ^ 0. This is in- 
deed consistent with what we observed; figure lTUI provides 
an example of this behavior. We confirm the values of 
AfsAQQ{nl) in Eq. (|46|l that can be immediately deduced 
from Tables Ill lIVI The analytic results for A{sAQo{nl) 
are thus found to be consistent with the numerical data 
for AfsGsE; the level of confirmation is 5-10 % [relative 
to AisAeoinl)] for P and D states (1 % for the 2P states, 
and 8D states not included, for the reason mentioned 
above), 3 % for F states, and 1% for G states. 

This represents an improvement over the accuracy 
of ^60(^1^^) obtained by the previous check. This im- 
provement comes evidently from the fact that the relative 
deviation of AfgGsE in Eq. lO from Afs^eo in Eq. 
is small over the whole range < Z < 110, compared to 
the relative deviation 



GsEjnlj, Za) 
Aeo{nlj) 



(48) 



of GsE [see Eq. ^] from Aeo{nlj) in Eq. © — with 
j = I + 1/2 01 j = I — 1/2. As a consequence, the un- 
certainty in the numerical evaluation of the limit of H47|) 
as ^ is relatively small. Figure [TUI shows an exam- 
ple of the smallness of the contributions to AfgGsE that 
go beyond AfgAgo- Moreover, we have observed that the 
higher the angular momentum I, the smaller the values of 
the deviation (|47|) , hence the stronger confirmation of our 
values of AfsAgo(n/) for high orbital angular momenta. 



VIII. SUMMARY OF RESULTS 

This paper contains results that are relevant to the 
self energy of a non-S electron bound to a point nucleus 
of charge number Z. We provided estimates and values 



(see also Ref. fo^) for the first two non-analytically-known 
contributions to the self-energy expansion (jSJ, namely 
the Bethe logarithm lnko{nl) and the so-called Aeoinlj) 
coefficient, which can be viewed as a relativistic Bethe 
logarithm. The main numerical results are contained 
in Tables HHIVl in Eq. ^ and Table IVnTl in Eq. ^ 
and in Eq. H37|) . We have also conjectured, in Sec. IVII 
that the relativistic Bethe logarithm AQQ{nlj) does not 
strongly depend on the principal quantum number n. In 
addition to this, we note that the orders of magnitude of 
AQQ{nli_i^2) and AQQ{n{l + 1)1+3/2) are the same (for a 
given set of quantum numbers n and I > 1), in Tables Q]- 
IIVI These results, taken together, yield in particular the 
best available approximations of the self energy in hy- 
drogen and light hydrogenlike ions, except for n = 1 and 
n — 2 levels ^l^il3] (see also Sec. IVII|) : such an approxi- 
mation can be obtained through Eqs. Q and 10. 

Calculating Ago has been a challenge since the seminal 
work of Bethe JJ on the dominant self-energy coefficients 
of S states [see Eqs. {Tjl and Q]. Details of the method 
we used were described in Sec. Illll and IIVI As discussed 
in Sec. IVIII including the coefficients Aqq reported in 
Tables H HrVl in a (truncated) expansion of the self energy 
improves its accuracy over a large range of nuclear charge 
numbers Z. 

We checked our calculations of Ago by both analytic 
and numerical means. The so-called e method, which 
we have employed (see Sec. IIIip . makes divergences ap- 
pear in the low- and high-energy contributions to Ago, 
as the scale-separating parameter e between these two 
contributions goes to zero. We have observed that, as re- 
quired, these divergences cancel when the two parts are 
added. Moreover, our calculations correctly reproduced 
the known lower-order coefficients A40 and Agi. We have 
also checked our results for Ago against numerical values 
of the self energy, and were able to confirm them by this 
independent method to the level of about 15 % (except 
for D3/2 states, as explained in Sect. lVIl|l . 

Obtaining results for Ago required extending (analyti- 
cally) the angular algebra developed for 2P states [13 to 
higher angular momenta. Techniques of numerical con- 
vergence acceleration of series 0, were instru- 
mental in evaluating the parts of Ago that could not be 
analytically calculated. The recent analytic calculations 
of Ref. 1231 enabled us to obtain with a high precision 
the self energy (0 of electrons with high {j > 3/2) angu- 
lar momentum, for various values of the nuclear charge 
number Z; the new calculations that we have performed 
required the use of massively parallel computers, and 
thousands of hours of computing time. (These numer- 
ical data, which have been used for the plots in Figs. 
ISHTUl will be presented in detail elsewhere [l^.) We 
have also collected the most recent available values of the 
self energy. This provided us with independent values of 
the Ago coefficients, extracted from the numerical self- 
energies, thus allowing us to check the analytic results 
presented in Tables Ultfvl f see SedVll. 

Severe cancellations appeared, between different con- 
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tributions to Aqq (in addition to the cancellation of the 
e-parameter divergences): for some of the atomic states 
investigated, the absolute magnitude of the Aqq coeffi- 
cients is as small as 10"'^, whereas the largest individual 
contribution to Ago, when following the classification of 
the corrections according to Refs. is of the order 

of 10~^ or larger for all atomic states discussed here (see 
also Tables and rvm. 

Future calculations of the Bethe logarithm lnko{nl) 
and of the relativistic Bethe logarithm AQo{nlj) could 
also fruitfully be compared to the estimates given by 
Eqs. 1(23), (EHl and (|27Il, and Table lvml The results pre- 
sented in this paper also allow one to perform checks of 
future exact self-energies obtained by numerical methods, 
by comparing their values to the three-term self-energy 
approximation (|42|l provided here for P and higher-/ 
states. The values of Aqq in Tables HI llVI can be of inter- 
est for analyzing the Lamb shift of highly-excited (high-n 
and high-/) electronic states in recent J,, JJJ and 

future high-precision spectroscopy experiments. The re- 
sults of Sect. IIVHVII also provide the best available self- 
energy approximation for many states nlj and nuclear 



charge numbers Z (see Sec. IVUfl : these approximations 
can for instance be useful in evaluating the contribution 
of QED effects in atoms 113, EH El or molecules [6^. 
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APPENDIX: LOCAL FITS 

This appendix describes a fitting procedure which is 
designed to extract "/oca/" numerical quantities from a 
set of data points, and to allow one to assess the numer- 
ical uncertainty associated to these quantities. A partial 
sketch of this procedure was first introduced in Ref. 65]. 
Here, "local" refer for instance to the evaluation of a per- 
turbation expansion about one abscissa; the purpose of 
the method presented here is to perform fits that are lo- 
cal to an abscissa of interest, as opposed to finding the 
best global fit of some data points. We thus used it in 
order to obtain asymptotic coefficients for AQQ{nlj) for P 
and D states in Sec. (see Table IVIIl|) . as well as the 
asymptotic expansion of the Bethe logarithm lnfco(ri/) in 
Eq. (|37|) — in these applications, the quantities evaluated 
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FIG. 11: This figure shows the lines going through a few pairs 
of successive data points 15011 — see also Fig. |5| Each of these 
lines is a local approximation to the curve underlying the data 
points. Each line yields an estimate of the limit 14911 of the 
data points as (this estimate is at the intersection of 

the line with the = axis). Fig. 1121 graphically displays 
these estimates. 
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FIG. 12: This figure shows the estimates of limit I49II ob- 
tained through the two-point fits of Fig. 1111 From this graph, 
we limit 1)49^ to be —0.0568(1), which is more precise than, 
and coherent with the value —0.057(1) obtained from the orig- 
inal data points H5()^ in Figs. |21 and 1111 The limit estimates 
are plotted along the vertical direction, while the abscissa as- 
sociated to an estimate is the average abscissa of the two data 
points of Fig. Illl that were used in producing it. 



are local to either n = oo or / = oo. This method can in 
principle be applied to many other problems that require 
local fits. 

In order to describe the local-fit procedure, we take the 
evaluation of the limit 



lim X lnA;o('^0 



(49) 



as an example — here, we have n = I + \ and \nkQ[nl) 
is the Bethe logarithm ©. This limit was evaluated as 
-0.056853(2) [see Fig. Hand Eq. (|S7|)]. 

Figures|2]andE]contain data points which are relevant 
to H49|l : we have plotted 



/^ X h\ko{nl) 



(50) 



as a function of (with values of the Bethe logarithm 
found in Ref. 0). The limit (|49|l can visually be esti- 
mated from the data points in Fig. |21to be —0.057(1). 

In order to improve over the estimate —0.057(1) 
for (|49|l . we fit (exactly) each pair of two consecutive 
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FIG. 13: From the lower to the higher curve: estimates of 
limit (1491 obtained through fits of the data points 1501 1 with 
polynomials of degree 1 (see also Fig. I12|l . 3 and 5 (see also 
Fig. [nil. Fitting the data points (|^ of Fig. (O with 1 to 6 
points yielded mutually coherent estimates of limit II49|I with 
an exponentially decreasing error. 
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FIG. 14: This figure shows estimates of limit 1491 obtained by 
fitting the data points 1501 in Fig. 1111 with fifth-degree poly- 
nomials (in The high relative stability of the estimates 
as — » allowed us to give the precise value —0.056853(2) 
in Eq. ^ for limit (gSJ. 



points H5U|I in Fig. [21 with a line, as depicted in Fig. 1111 
Each of the fitting hnes in Fig. ^2 gives an estimate of 
hniit (|49|l by extrapolation to = (intersection of 
the line with the = axis). Figure IT^ contains each 
of these estimates, as a function of the average abscissa 
of the two points that were used in obtaining it. Be- 
cause the curve in Fig. ^] is relatively flatter than the 
curve in Fig. 1111 we can estimate limit l|49|) with an im- 
proved uncertainty; thus, we deduce from Fig. E| the 
value —0.0568(1) for the limit H49I) that we are studying, 
which is coherent with the previous estimate —0.057(1). 

This better estimate -0.0568(1) of limit (gS)) can be 
further improved by continuing to increase the number p 
of data points (|50|l included in local fits of the data. Thus, 
for an increasing number p of data points, we fitted (ex- 
actly) each set of p successive points (|50|) in Fig. II II with 
a polynomial of degree p — 1 (linear combination of the 
functions 1, /^^,. . . , Z^^^^^^), and represented the value 
of the polynomial extrapolated tol^^ = as a function of 
the average abscissa of the p points. Fig. depicts this 
process. The plotted values are estimates of the limit (|49f) 
obtained with higher and higher-order (local) fits of the 
data points (|50|) . In Fig. ^1 the abscissa of each esti- 
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FIG. 15: This figure displays estimates of limit 1491 obtained 
by fitting the data points 1501 in Fig. Illl with a eigth-degree 
polynomials (in It should be compared to Fig. 1141 which 
gives a more accurate estimate of limit 1491 by fitting se- 
quences of only six data points. The accuracy of the local fits 
performed here first increases with the order of the local ap- 
proximations to the data points 1501 (see Fig. 1131 , and then 
eventually decreases (compare this plot to Fig. 1141 . 



mate is the average of the abscissas of the fitted data 
points (|50|l . We observed that the curves so obtained 
become exponentially flat, in the sense that their relative 
amplitude become exponentially smaller and smaller — 
until the uncertainties of individual estimates become 
important, as described below. This fact, which is illus- 
trated in Fig. ifT^ . allowed us to obtain more and more 
accurate estimates of limit (|49f) . 

The most accurate value that we obtained for 
limit H49() through the local-fit procedure described here 
is -0.056853(2) [see Eq. as is illustrated in Fig.lTl 

This limit was obtained by fitting each sequence of p = 6 
data points with a fifth-degree polynomial. Fits of the 
data points (|50|l with larger numbers of data points dis- 
play more irregular estimate curves; this can for instance 
be seen by comparing Fig. E| with Fig. 1151 

As we have seen above, the uncertainty in the fitted 
value can be evaluated by visually prolongating the fit- 
ting curves (i.e., curves such as those of Figs. I12II15|I . 
Another uncertainty must in general be taken into ac- 
count in order to obtain a reliable estimate for the fitted 
quantity: the uncertainty in the data points. All the 
curves presented in this appendix do contain error bars 
that reflect the uncertainties in the estimates of (|49|l that 
come from the uncertainties in the data points H50|l . We 
evaluated the uncertainty associated to each fit of p data 
points (|50|l by calculating three fits: a fit with the mid- 
dle values of the ordinates, a fit with the higher values, 
and a fit with the lower values; the three estimates of 
the fitted quantity (|49() obtained through this procedure 
define an estimate with an error bar (see, e.g.. Fig. I15|l . 
Other ways of estimating the uncertainty in the fit re- 
sult can be used; a good choice of uncertainty evaluation 
yields successive estimates of the fitted quantity that are 
compatible with a smooth curve of estimates [see, e.g.. 
Fig. El where the less precise estimates of limit (|49f) lie 
in the prolongation of the more precise values, which are 
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on the right of the plot]. 

One of the advantages of the local- fit method presented 
in this appendix is that data points that are located far 
from the abscissa of interest = 0, here) can fruitfully 
be used in evaluating the fitted quantity [limit H49|l . in our 
example]. Thus, as Fier. 1151 illustrates, data points H5()|l 
with "large" abscissas can yield more precise estimates 
of limit H49|l than data points with small abscissas. This 
behavior is particularly useful when data points in the 
region of interest have relatively large uncertainties. 

The procedure detailed in this Appendix also allows 
one to study the quality of lists of numerical results that 
should lie on a smooth curve, but whose coherence is not 



obvious through a simple inspection or plot of the values. 
In fact, curves such as those found in Figs. IT^HTKl can be 
very sensitive to small errors in a list of numerical values. 
We have not noticed such errors in the Aqq values of Ta- 
bles and ^1 while evaluating the asymptotic coefficients 
reported in Table IVlTll this provided an additional check 
of the values reported in these tables (see also Sec. IVIl) l. 

The local-fit method described here is not restricted 
to the asymptotic study of the Bethe logarithm that we 
have used as an example. In general, it can yield precise 
estimates of quantities that are local to a set of data point 
[such as limit (|49|l ] . including, for instance, perturbation 
coefficients of non-analytic expansions [e.g., Eq. (jsj]. 
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